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ABSTRACT 

The effect of five different C type grid geometries 
on the predicted heat transfer and aerodynamic perfor- 
mance of a turbine stator is examined. Predictions were 
obtained using two flow analysis codes. One was a finite 
difference analysis, and the other w*as a finite volume 
analysis. Differences among the grids in terms of heat 
transfer and overall performance were small. The most 
significant difference among the five grids occurred in 
the prediction of pitchwise variation in total pressure. 
There was consistency between results obtained with 
each of the flow analysis codes when the same grid w*as 
used. A grid generating procedure in which the viscous 
grid is embedded within an inviscid type grid resulted 
in the best overall performance. 
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INTRODUCTION 

One of the goals for Navier-Stokes turbomachinery 
analysis is the accurate predictions of blade row losses 
and turbine blade heat transfer. It has been shown by 
a number of investigators that the characteristics of the 
grids used can significantly affect the predicted results. 
Davis, Hobbs and Weingold(1988) investigated the ef- 
fects of blade-to-blade grids of different densities on the 
prediction of compressor performance. Recently, atten- 
tion has also been given to the desirability of using grids 
orthogonal to the blade surfaces. Many applications us- 
ing C-type grids have used grids generated in a manner 
similar to that proposed by Sorenson (1980) in which 
the points on the cut line have the same coordinates 
for both the upper and lower portions of the grid. It 
was shown by Arnone, Liou, and Povineili(1992) that 
C type grids can be highly skewed, especially for highly 
turned modern turbine stator blades. The skewing is 
most evidenced when a matching condition is enforced 
on the cut line from the trailing edge to the downstream 
boundary. If the matching condition is removed the re- 
sult is less skewed grids. The cut line divides the C 
grid, and extends from the trailing edge to the down- 
stream boundary. Arnone et al.(1992) proposed the use 
of C-type grids for turbine applications in which the re- 
quirement for common coordinates along the cut line 
is not enforced. This procedure allows for greater grid 
orthogonality at the blade surface, since the number 
of grid cells on the upper and lower surfaces are not 
necessarily equal. However, this procedure requires in- 
terpolation along the cut line, rather than averaging 
the solution at a common point. Micklow, Shivaraman, 
and Li( 1993) took this process further, and presented 
results where the requirement for matching points was 
not enforced along the outer boundary of the blade-to- 
blade grid, as well as along the cut line. 

Yeuan, Hamed, and Tabakoff(1993) used non- 
periodic H grid to analyze the performance of a tur- 
bine cascade. Turner, Liang, Beauchamp, and Jen- 
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nions(1993) advance arguments for the use of modified 
H grids, called I grids, to improve flow field calculation 
accuracy. The H grids were modified so as to be more 
orthogonal to the blade surface, and were recommended 
for use in calculating the performance of highly turned 
turbine blades. Because of the rapid changes in heat 
transfer in the leading edge region C or 0 type grids 
are advantageous in comparison to H or I type grids 
for use in Navier-Stokes heat transfer analyses. C type 
grids have the advantage relative to 0 type grids in that 
only a single grid is needed for a blade row passage. 0 
type grids are generally embedded within a H type grid 
for cascade analyses. The work that is reported in the 
present paper is concerned with determining the effects 
of different C type grids on turbine blade heat transfer, 
as well as turbine blade aerodynamic performance. 

The issue of grid size is especially important w hen 
one considers that the goad of code development work 
is to achieve accurate three-dimensional Navier-Stokes 
solutions while utilizing a reasonable amount of com- 
puter resources. It has been shown by Boyle and 
Giel(1992) that over 50 spanwise grid planes are nec- 
essary to achieve grid independent heat transfer results 
for a typical turbine blade. If one wishes to limit the size 
of routine three dimensional Navier-Stokes calculations 
to around a half million grid points, then blade-to-blade 
grids would have a maximum of 10,000 points. 

The choice of the best blade-to-blade grid is facil- 
itated by comparisons with experimental data. Arts, 
Lambert de Rouvroit, and Rutherford(1990) presented 
data for the midspan region of a turbine stator with an 
exit flow angle of 7 4°. Heat transfer, pitchwise variation 
in total pressure, as well as overall stator performance 
data were given. These data were used by Harasagama, 
Tarada, Baumann, Crawford, and Neelakantan(1993) 
to compare different approaches for the prediction of 
turbine blade heat transfer using boundary layer meth- 
ods. Luo and Lakshminarayana(1993) used the same 
experimental data to show the validity of their method 
for predicting the flowfield and heat transfer in a tur- 
bine passage using a Navier Stokes analysis. 

In the present work blade surface heat transfer as 
well as aerodynamic performance are examined for five 
different baseline grids. Additional grids were gener- 
ated to examine other possible grid effects. Three of 
the five baseline grids have a matching condition im- 
posed along the cut line. Among these three grids two 
are similar, and differ only in the procedure used to 
generate them. One was generated using a single grid 
solution generated using the GRAPE code of Soren- 
son(1980). The other was generated using the code of 
Arnone et al.(1992) in which the viscous grid is gen- 


erated by embedding a near wall grid within a sparse 
inviscid grid. The third grid that is periodic along the 
cut line differs from the other two in that the grid in the 
wake region expands to have nearly uniform pitchwise 
spacing at the downstream boundary. The two grids 
without the matching condition differ from each other 
downstream of the blade row. One extends the periodic 
boundary in a straight line, so that there is little flow 
across the periodic boundary downstream of the blade. 
The other curves the grid, so that there are large flows 
across the periodic boundary downstream of the blade. 

The choice of grid could be significantly affected by 
the Navier-Stokes solution methodology. The Navier- 
Stokes analysis were done using a finite difference code 
developed by Chima(1987), and a finite volume code de- 
veloped by Arnone et al.(1992). Since the choice of cal- 
culation procedure is affected by far more factors than 
are explored in this paper, the emphasis is on deter- 
mining the best grid approach for each code. Similarly, 
since the emphasis is on grid effects, a simple approach 
is taken to turbulence modeling. A mixing length tur- 
bulence model is used. 


DESCRIPTION OF ANALYSIS 


Grid characteristics. All five grids were of of the same 
size, 177 x 53. This size was chosen based on re- 
quirements for three-dimensional Navier-Stokes analy- 
ses. Certainly, two-dimensional Navier-Stokes results 
could be obtained in a reasonable CPU time for grids 
of larger sizes. However, based on the required num- 
ber of spanwise grids for three dimensional heat trans- 
fer analyses, the CPU time and core memory required 
for three-dimensional analyses might be excessive for 
blade-to-blade grids significantly greater than 177 x 53. 
All grids, except grid A, had the inlet boundary 0.55c r 
ahead of the vane. Grid A began 0.75c x in front of 
the vane. All grids had their downstream boundary at 
0.84c r behind the vane trailing edge. The downstream 
boundary w*as chosen so that the measurement sta- 
tion used by Arts, Lambert de Rouvroit, and Ruther- 
ford(1990) was midway between the trailing edge and 
the downstream boundary. The near wall spacing was 
held constant at 4 x 10“ 5 c r for all baseline grids. The 
effect of variations in the near wall spacing on predicted 
results is examined. 


2 



MATCHING & NON-MATCHING WAKE CUT LINES 




BEHIND VANE 




Fig. 1 - Description of grids used in analysis 
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The five grids used in the analysis are illustrated 
in figure 1, and will be referred to with the labels A 
through E. Table I give the number of grid spacing 
used in each of the baseline C-type grids. Grid A was 
generated using the GRAPE code of Sorenson(1980). 
This code obtains a solution to a two-dimensional Pois- 
son partial differential equation in order to generate a 
C-type grid. Grid B was generated using the code de- 
veloped by Arnone et al.(1992), and as can be seen in 
figure 1, this grid is very similar to grid A. The princi- 
pal difference is that a two step grid generation process 
was used to generate grid B. First a coarse grid, such as 
might be used in an inviscid analysis is generated. This 
initial grid has relatively few lines in the blade-to-blade 
direction. The grid used for the viscous calculations is 
obtained by embedding a fine grid over a few cells in 
the near wall region. The remaining three grids were 
generated using the same procedure as grid B. Grid C 
differed from grid B in that the grid lines for small val- 
ues of rj do not maintain a constant spacing from the 
cut line between the trailing edge and the downstream 
boundary. For grid B these grid lines are maintained 
at & uniform spacing between the trailing edge and the 
downstream boundary. Grids D and E differed from 
grid C in that these two grids were not periodic along 
the cut line between the trailing edge and the down- 
stream boundary. As can be seen in figure 1, grids that 
are not periodic along the cut line can more easily be 
made orthogonal to the blade surface. One potential 
disadvantage of this approach however, is that the flow 
solution has to be interpolated along the cut line. Grid 
E differs from grid D in that downstream of the blade 
the periodic boundary is curved so that at the down- 
stream boundary the tj = constant grid lines are parallel 
to the axial direction. Grids A through D are extended 
downstream of the blade at a constant angle, which is 
dose to the trailing edge angle. When these grids are 
used, there is only a single wake in the flow field. How- 
ever, when grid E is used, multiple wakes are present in 
the flow field. The number of wakes depend on the flow 
angle and the distance between the trailing edge and 
the downstream boundary. At any given axial location 
there is only a single passage flow field. Grids Dn and 
Ea in table I are similar to grids D and E respectively, 
but with more circumferential grid lines. 

Table I. - Distribution of circumferential grid. 



Grid 

No. of increments 

A 

B 

C 

D 

E 

Dn 

Ea 

Wake-pressure side 

40 

40 

40 

56 

56 

88 

128 

Blade- pressure side 

48 

48 

48 

32 

32 

32 

32 

Blade-suction side 

48 

48 

48 

64 

64 

96 

96 

Wake-suction side 

40 

40 

40 

24 

24 

24 

64 


Flow analysis. Two flow analysis codes were run for 
each of the five grids. One analysis used a finite dif- 
ference approach, and was developed by Chima(1987). 
In the discussion which follows, results obtained using 
this code are labeled as flow code FC. The other anal- 
ysis, developed be Arnone et al.(1992), used a finite 
volume approach, and employed a multi-grid solution 
scheme. Results obtained using this code are labeled as 
flow code FA. While the discretization is different be- 
tween the two analyses, both used a time marching ap- 
proach with an explicit four stage Runge-Kutta scheme 
to solve the differential equations. Both also employed 
implicit residual smoothing. 

The turbulence model used in both flow analyses 
is a variation of the model developed by Baldwin and 
Lomax(1978). The main difference between the model 
used and the Bald win- Lomax model is in the prediction 
of the transition location. The transition model given 
by Mayle(1991) was used. In this model the location 
of the start of transition as well as the length of tran- 
sition are given as functions of Reynolds number and 
turbulence intensity. It was found that estimating the 
local freestream turbulence intensity rather than just 
using the inlet freestream turbulence intensity gave a 
better estimate of the vane heat transfer for the cases 
examined. The local freestream turbulence intensity 
was estimated assuming that the velocity fluctuations 
remain constant through the the passage. This assump- 
tion results in: 

Tus = TuisletUinlbt/Us 

The local velocity, Us 9 was calculated from the isen- 
tropic relationship and the local static- to- inlet total 
pressure ratio. 

RESULTS 

Heat Transfer Comparisons. Figures 2 and 3 show the 
variation in heat transfer among the different grids for 
the two flow solution codes. Also shown in each figure 
are the experimental data of Arts et al.(1990). These 
comparisons are for Rt* = 1.16 x 10 6 , Tvinlet = 1%, 
and M-x — 0.84. Because of the low inlet turbulence, 
transition did not occur even though the Reynlods num- 
ber is fairly high. The Stanton number predictions 
shown in figure 2 show very high heat transfer just near 
the trailing edge. These results were obtained by forc- 
ing the flow turbulent at the tangent point of the vane 
trailing edge circle. Since the focus of this work is on 
the effects of different grids on predicted results, it is 
useful to show comparisons for a laminar flow case. For 
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St x 1000 



Fig. 2 Stanton number comparisons for Re 2 = 1.16 x 
10 6 ,Tu = 1%, \f 2 = 0.84, Flow code FA. 



Fig. 3 Stanton number comparisons for Re 2 = 1.16 x 
10 6 ,Tu = 1%, M 2 = 0.84, Flow code FC. 



Surface distance, s/c 

Fig. 4 Effect of transition model assumptions on pre- 
dicted heat transfer, Re 2 — 2.1 x 10 6 ,Xu = 6%, M 2 — 
0.92, Flow code FC, grid C. 

laminar flow, questions regarding the implementation of 
the turbulence model do not arise. For each of the flow 
analyses, the variation in blade heat transfer among the 
five grids is small. Exept very close to the trailing edge, 
both codes give nearly the same level of heat transfer for 
the same grid. In flow code FA the boundary layer was 
set to be fully turbulent at the beginning of the trailing 
edg circle. The agreement with the experimental data is 
generally good, and for this case there is little evidence 
to prefer one grid over another. 

The heat transfer predictions are in good agreement 
with the experimental data for the low turbulence in- 
tensity case. Comparisons will next be made for a case 
with high turbulence intensity, and high Reynolds num- 
ber. The results of different approaches to modeling the 
effects of a high turbulence intensity will be examined 
prior to showing the grid effects. Figure 4 shows the 
effects of transition model assumptions on vane surface 
heat transfer for the highest Reynolds number, 2.1 x 10 6 , 
and highest turbulence intensity, 6%, tested by Arts et 
al(1990). The exit Mach number was 0.92. Results 
are shown for a single grid, (grid C), and flow solution 
code, (FC). The purpose of this figure is to illustrate the 
significance of the transition model for the heat trans- 
fer predictions. The model for the start of transition 
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proposed by Mayie(1991) more accurately predicts the 
start of transition when the turbulence intensity is ad- 
justed to account for the local inviscid velocity. This 
is especially evident on the suction surface. Comparing 
the slope of the predicted heat transfer with the ex- 
perimental data shows that the length of transition on 
the suction surface is not well predicted using the in- 
termittency model proposed by Mayle(1991). The pre- 
dicted transition length is shorter than that indicated 
by the experimental data. Heat transfer predictions 
are shown for two additional transition length models. 
These are the transition length models of Simon and 
Step hens( 1991) and Simon(1994). The model of Simon 
and Stephens was developed for zero pressure gradi- 
ent flows, while Simon’s model was developed for flow's 
with favorable pressure gradient. The use of Simon’s 
transition length model results in good agreement with 
the experimental data. Even though suction surface 
transition began close to the uncovered portion of the 
suction surface, there w'as a favorable pressure gradient 
at the start of transition. These results illustrate the 
sensitivity of the heat transfer results to transition as- 
sumptions. Because of the high freestream turbulence 
intensity, pressure surface transition occurs close to the 
leading edge. The heat transfer on the pressure surface 
is accurately predicted. 

Figure 5 shows the effect of different model assump- 
tions for predicting the effect of freestream turbulence 
on blade surface heat transfer. Heat transfer predic- 
tions are shown using the model of Forrest (1970), as 
w'ell as when no augmentation due to freestream turbu- 
lence is assumed. This model for calculating an eddy 
viscosity due to freestream turbulence was applied only 
prior to transition. It was applied in an analogous fash- 
ion to the Baldwin- Lomax eddy viscosity model, which 
is used after transition occurs. The Baldwin-Lomax 
model for turbulent eddy viscosity is a two layer model. 
Predictions are show'n w’hen the augmentation model is 
applied only in the inner layer; when it is applied to 
both layers; and w'hen it is held constant in the outer 
layer. This model for the augmentation of eddy viscos- 
ity is utilized only prior to transition. The best over- 
all agreement with the experimental data is achieved 
w'hen the turbulent eddy viscosity in the outer layer is 
held constant at the inner layer value determined at the 
point where the two layers meet. 

A comparison of the heat transfer predictions for 
the five different grids is given in figures 6 and 7. Since 
this is a relative comparison, only Mayle’s transition 
model is utilized. The use of this model results in a rel- 
atively long portion of the suction surface being turbu- 
lent, and allows for the comparisons for different grids 


to be done for a case with turbulent flow on both pres- 
sure and suction surfaces. In this case the pressure 
surface is mostly turbulent, and half of the suction sur- 
face is turbulent. The variation in blade surface heat 
transfer among the five grids is small. The results for 
this test case are similar to those shown in figures 2 and 
3 for a case where the vane surface heat transfer was 
almost entirely laminar. 

The variation in predicted blade surface heat trans- 
fer between the tw'o flow' codes for a given grid is some- 
w'hat larger than the variation due to different grids 
for either of the two flow analyses. Part of the differ- 
ences can be attributed to differences in modeling the 
turbulent eddy viscosity. In the leading edge region, 
and prior to transition, flow code FC used the model 
of Forrest (1977), to augment the laminar viscosity for 
the effects of freestream turbulence. This model w'as 
applied to both the inner and outer regions, since this 
assumption would magnify any differences due to dif- 
ferent grids on the heat transfer. Flow code FA does 
not increase the viscosity to account for the effect of 
freestream turbulence prior to transition. Both flow 
analyses overpredict the suction surface heat transfer 
when the flow is fully turbulent. 



Surface distance, s/c 

Fig. 5 Effect of model assumptions for freestream 
turbulence on predicted heat transfer, Re 2 = 2.1 x 
10 6 , Tu = 6%, M 2 = 0.92, Flow code FC, grid C. 
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Fig. 6 Stanton number comparisons for Re 2 = 2.1 x 
10 6 ,Tu = 6%, Mo = 0.92, Flow code FA. 



Surface distance, s/c 


Fig. 7 Stanton number comparisons for Reo = 2.1 x 
10 6 ,Tu = 6%, Mo = 0.92, Flow code FC. 


The results shown in figures 2 and 3 are for Re 2 = 
1.16 x 10 6 , while those in figures 6 and 7 are for a 
Reynolds number nearly twice as large. Since the same 
grids were used, the value of y* for the results in fig- 
ures 6 and 7 is greater. The maximum value of y* for 
the high Reynolds number case was 3.8. The use of the 
largest value of yj 1 " which gives accurate results mini- 
mizes the amount of near-wall grid stretching, and also 
allows for the smallest size grid in the blade-to-blade 
direction. The effect of reducing the near wall spacing 
is shown in figure 8 for grid D, and flow code FC. There 
was virtually no change in the heat transfer as the near 
wall spacing was reduced by a factor of 4. Consequently, 
the relatively large near wall spacing used for the high 
Reynolds number analysis is not a cause of the differ- 
ences between the measured and predicted vane heat 
transfer. 



Surface distance, s/c 


Fig. 8 Grid spacing effects on heat transfer. Re 2 = 
2.1 x 10 6 , Tu = 6%, M 2 = 0.92, Flow code FC, grid C. 

Reducing the near wall spacing did not result in 
better agreement with the experimental heat transfer 
in the fully turbulent region. To determine if this dis- 
agreement was related to the grid properties, the effect 
of improving the orthogonality of the grid lines to the 
vane surface was examined. Figure 9 shows the angle 
that the grid lines make with the vane surface. If the 
lines were purely orthogonal, the angle would be ninety 
degrees. Grids B and C make identical angles with the 


blade surface, since they differ only downstream of the 
blade. The embedded grids, (B through E), have grid 
lines which maintain a constant angle for the region be- 
tween the blade surface and the first circumferential in- 
viscid grid line. For the viscous embedded grid this an- 
gle is maintained over twenty near wall spacings. Grid 
A is solved as a single grid, and the angle of the grid line 
with the blade surface continuously varies away from 
the wall. While grids B and C have the highest depar- 
ture from orthogonality, it should be noted that they 
are both periodic in the wake. Enforcing grid matching 
along the cut line results in embedded grids w’hich have 
a higher degree of nonorthogonality at the blade sur- 
face than either grids D or E. W r hile grids that match 
along the cut line can be generated using the procedure 
discussed by Arnone et al.(1992), they did not recom- 
mend doing so. Figure 10 show's the angle of the grid 
lines with the blade surface for grids D and Dn. As 
seen in table I, grid Dn is significantly larger in size 
than grid D. The use of a larger grid allows for greater 
orthogonality at the blade surface. However, figure 11 
shows that the more orthogonal grid did not result in a 
significantly different heat transfer distribution. 



Fig. 10 Angle between grid line and blade surface for 
grids D and Dn. 



Surface distance, s/c 


Fig. 9 Angle between grid line and vane surface. 



Surface distance, s/c 


Fig. 11 Stanton number comparisons - grids D and Dn, 
Re 2 = 2.1 x 10 6 , Tu = 6%, M 2 = 0.92, Flow code FC. 
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Table II. - Overall loss coefficient. 


Total pressure distribution. Arts et al.(1990) mea- 
sured the pitch wise variation in total pressure at 0.42 
axial chords behind the vane. Figures 12 and 13 com- 
pare the pitchwise variation in total pressure for the five 
grids and the two flow solvers with the experimental 
data. For clarity of presentation, the minimum pres- 
sure was taken as the abscissa origin. Therefore, any 
variation in flow angle among the grids is not reflected 
in these figures. The relative grid effects are similar for 
both flow solvers. For grids A through D flow solver FA 
resulted in minimum total pressures significantly less 
than flow solver FC. The results obtained using grids A 
through D showed the predicted minimum total pres- 
sure lower than the experimental data. The width of 
the predicted wake is also smaller than the experimen- 
tal wake. This indicates that the amount of diffusion 
of the wake predicted by the turbulence model is too 
small. Results obtained with flow solver FC and grid 
A showed a lower total pressure in the freestream re- 
gion than with grids B, C, or D. In this region the data 
show no loss in total pressure. Results obtained with 
flow solver FA and grids A through D showed no loss 
in tot ad pressure in the freestream region. 

Grid E is bent so that the grid is aligned with the 
axis at the downstream boundary. This bending results 
in a relatively coarse grid on the suction side of the cut 
line. The reasonably good agreement between the pre- 
dicted and measured wakes for grid E may have been 
fortuitous, The numerical diffusion resulting from the 
coarse grid may have helped to give nearly the correct 
total diffusion. To verify this hypothesis, an additional 
bent grid, (Ea), was generated. Table I shows that the 
number of points in the wake was increased consider- 
ably. Figure 14 compares the wake profiles for the two 
grids. The narrowing of the wake, and lowering of the 
minimum total pressure, in going to the finer grid is 
evidenced. 

A comparison of the wake profile using a two- 
equation turbulent kinetic energy turbulence model 
with the wake profile using the Baldwin-Lomax tur- 
bulence model is given in figure 15. The results were 
obtained using grid D, and flow code FA. The two- 
equation model is that of Chien(1982), and details of 
its implementation into the flow code FA are given by 
Ameri and Arnone(1992). The depth of the wake using 
Chien’s model is nearly the same as the wake depth cal- 
culated using the Baldwin-Lomax model. The width of 
the w’ake is somewhat wider using Chien’s model. The 
wider wake region using Chien’s model is due to the 
the model giving transition on both surfaces close to 
the leading edge. As shown in figure 2, the Baldwin- 
Lomax model did not result in an early transition. 



Grid 


A 

B 

c 

D 

E 

Dn 

Ea 

Flow code 

e 

FA 

.039 

.029 

.030 

.030 

.054 

.028 

.032 

FC 

.035 

.026 

.023 

.029 

.024 

.027 

.024 



Experimental - Arts et al. 

0.029 


Overall lo6s. Table II shows the loss coefficient for 
the different grids obtained with the two flow solvers at 
an exit isentropic Mach number of 0.85. Except for the 
results obtained with flow code FA and grid E, the high- 
est calculated loss was achieved using grid A. The losses 
for grid E are different between flow codes FA and FC. 
But, based on the previous discussion, the results us- 
ing this grid are not expected to be accurate. For both 
flow codes the loss levels for the other grids are rea- 
sonably close to the experimental measurements. The 
losses calculated using flow code FC are slightly lower 
than those calculated using flow code FA. Part of this 
difference in loss is due to the blade surface boundary 
layers being tripped near the trailing edge when flow 
code FA was used. The results for grid D are in best 
agreement with the experimental data. 

The overall loss increased from 0.030 to 0.057 when 
Chiefs model was used. Part of the increase was the 
result of the boundary layers being turbulent. A fully 
turbulent calculation using the Badd win- Lomax model 
gave a loss coefficient of 0.041. 



Fig. 12 Pitchwise variation in total pressure, Re 2 = 
1. x 10 6 ,Tu = 1%, M 2 = 0.85, Flow code FA. 
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Pitchwise distance 


Fig. 13 Pitchwise variation in total pressure, 
Re 2 = 1. x 10 6 ,Tu=l%, A/2=0.8 o, Flow code FC. 



Pitchwise distance 

Fig. 14 Grid density effect on total pressure distribu- 
tion, Re 2 = 1. x 10^, Tu = 1%, A/o = 0.85, Flow code 
FC, 

CONCLUDING REMARKS 

The results of this investigation showed that the 
principle effect of different grid geometries examined 
was in the pressure distribution behind the vane. Four 
of the baseline grids resulted in an excessive decrease in 
total pressure at the center of the wake. The calculated 



Fig. 15 Wake profiles for Baldwin-Lomax and Chien 
turbulence models. Grid D, Flow code FA. 

wake was deeper, but more narrow than the experimen- 
tal data. The implication of this is that the turbulence 
model used gave insufficient physical diffusion in the 
wake region. The fifth baseline grid, which was rela- 
tively coarse in the wake region, gave good agreement 
with the experimental data for the minimum pitchwise 
total pressure. The numerical diffusion caused by the 
coarse grid resulted in a more nearly correct wake pro- 
file for this grid. A similar, but denser grid(Ea), which 
resulted in less numerical diffusion, gave wake profiles 
similar to the other four baseline grids. The relative 
effect of different grids was the same for either of the 
two flow codes used. 

The effect of different grid geometries on the vane 
surface heat transfer was small. For the low turbulence 
intensity case the agreement with the experimental data 
was good. For the high turbulence-high Reynolds num- 
ber case the agreement was influenced by the choice 
of model to determine transition length and intermit- 
tency. The start of transition was better predicted by 
Mayle’s transition model when the local turbulence in- 
tensity was taken as a function of the blade pressure dis- 
tribution. Simon’s transition length mode resulted in 
good agreement with the data. The choice of model to 
account for freestream turbulence significantly affected 
the predicted heat transfer. Forrest’s model gave rea- 
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sonably good agreement with the data when the aug- 
mented eddy viscosity was limited in the outer region 
of the two-layer Baldwin-Lomax turbulence model. 

The overall loss distribution was not strongly af- 
fected by the choice of grid geometries. Using either 
flow analysis the overall loss was reasonably well pre- 
dicted. 

In terms of the overall results grid D appears to be 
the best choice. This grid has a non-matching condition 
along the cut line, and the cut line is extended in a 
straight line. If it is desired to maintain a matching 
condition along the cut line either grids B or C appear 
acceptable. Results obtained using either of these grids 
were nearly identical. From the standpoint of speed 
of convergence, there w r as no reason to prefer one grid 
over another. The number of time steps to obtain a 
converged solution was about the same for all of the 
grids. Surface pressure distributions were also similar 
for each of the baseline grids. 
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